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Abstract. We study the dynamics of the Rabi Hamiltonian in the medium coupling 
regime with \g/u>\ ~ 0.07, where g is atom-field coupling constant, uj is the field 
frequency, for the quantum state with average photon number h ~ 10 4 . We 
map the original Hamiltonian to an effective one, which describes a tight-binding 
chain subjected to a staggered linear potential. It is shown that the photon 
probability distribution of a Gaussian-type state exhibits the amplitude modulated 
Bloch oscillations (BOs), which is a superposition of two conventional BOs with a 
half-BO-period delay between them and is essentially another type of Bloch-Zener 
oscillation. The probability transition between the two BOs can be controlled and 
suppressed by the ratio gy/K/co, as well as in-phase resonant oscillating atomic 
frequency fi (t), leading to multiple zero-transition points. 
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1. Introduction 

The interaction of matter and light is one of the fundamental processes occurring in 
nature. Its elemental constituent is the coupling between single atoms and photons. The 
quantum Rabi model considers a two-level atom coupled to a quantized field, describing 
the simplest interaction between quantum light and matter. It was used to describe 
the interaction between a rapidly varying, weak magnetic field and nuclear spin BIZI. 

Now it applies to a variety of physical systems with the two-level atom, being trapped 
ions, quantum dots, and superconducting qubits. Such systems can be exploited as a 
building block for quantum information processing and other potential applications to 
future quantum technologies [3j. 

Although the quantum Rabi model seems relatively simple and was declared to be 
solved recently M, the dynamics are actually quite complicated, depending on the 
system parameters and the initial state [6|. The rotating wave approximation (RWA) 
is justified in the strong coupling regime (but the average photon number can not be 
too large) and results in the Jaynes-Cummings (JC) model [7], which predicts non- 
classical phenomena, such as revivals of the initial excited state of the atom PEI Ho]. 

It is well known that an atom-cavity system undergoes the Rabi oscillations. In the 
case of the single cavity mode being initially prepared in a coherent state with large 
average photon number, although the oscillations experience the collapse-and-revival, 
the photon probability distribution remains unchanged in the framework of the JC 
model. The question is whether it is true if the original Rabi model is considered, or 
what will be observed in experiment. With strong couplings in the solid-state-cavity- 
QED m, the failure of RWA has relighted the interest on the quantum Rabi model 

mmmmmmmmm- 

In this paper, we study the dynamics of the Rabi Hamiltonian in the medium 
coupling regime with \g/co\ ~ 0.07, where g is atom-field coupling constant, to is the 
field frequency, for the quantum state with average photon number n rs-/ 10 4 . The 
effective coupling constant becomes \gVn\ ~ 7cn, which goes beyond the RWA. Within 
this parameter regime, the original Hamiltonian is mapped to an effective one, which 
describes a uniform tight-binding chain subjected to a staggered linear potential. An 
approximate solution suggests that the photon number distribution of a Gaussian- 
type state exhibits the amplitude modulated Bloch oscillations (BOs), which is a 
superposition of two conventional BOs with a half-BO-period delay between them. This 
phenomenon essentially belongs to Bloch-Zener oscillation in the strong tunneling limit 
20- EH [22] 123,' EH EH]. We find that the probability transition between the two BOs can 
be controlled and suppressed by the ratio g^/h/co, as well as resonant oscillating atomic 
frequency leading to multiple zero-transition points. A numerical simulation of 

dynamics in the Rabi model confirms this prediction. 

This paper is organized as follows. In Section [2] we propose the effective 
Hamiltonian in the concerned parameter regime. In Section [3] we present the 
approximate solution of the effective Hamiltonian, based on which the dynamics of 
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a local state is investigated. Section [I] is dedicated to a numerical simulation in the 
original Rabi model, focusing the control of the probability transition between two BOs 
by the atomic transition frequency. Finally, we give a summary and discussion in Section 

13 


2. Model and equivalent Hamiltonian 


We start by considering the single-mode atom-cavity model whose Hamiltonian can be 
written as 

H = ud)a + ^a x + gcr z {a} + a), (1) 

where oj and hi are the held and atomic transition frequencies, respectively, and g is 
the coupling constant, af (a) is the creation (annihilation) operator of the light held, 
while & x = |e) (e| — | g) (g |, d z = \e) (g\ + | g) (e| are atomic operators, where | g) and \e) 
denote the ground and excited atomic states, respectively. There is an approximation 
which has been developed, the rotating-wave approximation (RWA) [26, [27j under the 
condition |(?\/h| -C u. 

In this paper, we consider the case with the medium coupling regime with \g/oo\ ~ 
0.07 and the average photon number n ~ 10 4 . Such a parameter regime is accessible 
in experiments [28) [29, [30]. The aim of this paper is to investigate the dynamics of a 
particular initial state with a Gaussian-type photon number distribution, which allows 
us to obtain the approximate analytical result from the Hamiltonian ([TJ). 

Since Hamiltionian (JTJ) is parity-conserving for the excitation number 

A f = a^a + — (a x + 1), (2) 

it can be written in two independent equivalent Hamiltonian (Hf q ) in the bases 
with even (odd) excitation number: 


Heq — 9 ^ ^ 1 (N)a (* + 1| + H.C.) 

i =0 



(_f)7A+* 

2 


G + iu 




( 3 ) 


i=0 L J 

where A = e, o and q e = 1, q 0 = 0. This is schematically illustrated in Fig. [D 
The equivalence between H and Hf q is based on the mapping of the corresponding 
basis: for even invariant subspace we have {|g, 0), |e, 1), |g, 2), |e, 3), |g, 4) , - - •} 
—{|*) e }> where i = 0,1,2, •••, while for odd invariant subspace we have 

{|e, 0), |g, 1), |e, 2), |g, 3), |e, 4), • • •} —» {|i) 0 } |6] . Note that Hf q is a standard tight- 
binding chain with coordinate-dependent nearest neighbor (NN) hopping strength and 
on-site potential. It is easy to find out that, within the region of large excitation- 
number limit i » 1, when the initial state is local near n ~ 10 4 , we can take truncated 
approximation to Hf q . In this case, the expansion of the coupling strength in equation 
(j3J) is \fi = \/h + — g,r "/2 +.... We take the first term only, then y/i y/fi: therefore, 
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Figure 1 . (Color online) The atom-cavity level structure used for the mapping to a 
tight-binding chain, and the implementation of Bloch oscillations in the photon-number 
space. Atom-cavity level diagram showing the lower-energy states for a two-level atom 
of transition frequency fl coupled (with single-photon Rabi frequency g) to a single 
mode cavity field of frequency cu. with f2 = w. (a) and (b) correspond to even- and 
odd-parity cases, respectively. The transition strengths are in units of g. (c) For large 
photon region n 1, both even- and odd-parity cases can be mapped to a tight- 
binding chain with uniform NN hopping strength (y/n + 1 ~ y/n) and linear potential 
(as represented by gray shadow) approximately, which allows the occurrence of Bloch 
oscillations. 


H, e A q can be approximately equivalent to a tight-binding chain with uniform NN hopping 
strength and staggered linear potential. In the following, we neglect the label A in the 
effective Hamiltonian for simplicity and demonstrate that such a discrete system admits 
the existence of a robust BO [31] 32 133] . 

3. Approximate formalism 

In this section, we will introduce an approximate formalism for the approximate solution 
of the effective Hamiltonian and investigate the dynamical behavior as an application. 

3.1. Effective Hamiltonian 

We rewrite the effective Hamiltonian in the form 
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o N 

Hi ( 6 ) 

i =0 

where n is the average number of photon and can be regarded as a constant in the 
context of the problem we concern in this paper. This Hamiltonian is equivalent to (J3]) 
for the dynamics of a local state within the region ? > 1. Notice that the hamiltonian 
H e s is nothing but the tight-binding Hamiltonian to describe a single particle subjected 
to a staggered linear potential, which has been well studied in previous literatures [20]. 

It is demonstrated that the Bloch-Zener oscillations with a single period occurs in the 
parameter regime \g/oj\ ~ 0.2 and 12 = 6.734cu. In this paper, we revisit the same 
model within the parameter regime \g/u)\ ~ 0.07 and Q = oj, which is accessible 
in experiment |28j j29[ 30]• We will show that the dynamics undergoes Bloch-Zener 
oscillations characterized by two periods. 

In the absence of term Hi, it is a standard model which admits the existence 
of Wannier-Stark localization and Bloch oscillations, and has been studied extensively 
[33, [35| [36,, [37]. Now the question is: what happens to the dynamics of the model in 
Eq. ([4]). The basic idea is to consider the term H\ as a perturbation. According to the 
theory of Wannier-Stark localization [3T1 1 32] , the solution of Hq has the form 


H 0 | fj m ) = E° m \if m ), 

(7) 

o N 

E m = um - —uj, 

(8) 


and 


I '4 > m) 

I k) 


1 

Vn 

i 

Vn 


ex p 

k 

N 



mk H-sin k 

2 


J2e ikl \l),m = 0,±l,±2,.. 

i=i 


I k) 


In real space, it reads 

27V 

ItOn) ^ [ J l—m 
1=1 


io. 


(9) 

( 10 ) 


( 11 ) 


where L = —\g\ffijijj , is the spatial extent of a single Bloch wave packet oscillation [32], 
and g is negative. 

In order to examine the effect of Hi on the dynamics of the system, it is convenient 
to work in the interaction picture. The propagator of the whole Hamiltonian H in the 
interaction picture is obtained by the unitary transformation 


H (■ t) = e iHot Hie- iHot . 


( 12 ) 


The propagator is 

U 1 (t, 0) = Texp 



H t (; t') dt' 


) 


(13) 
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where T is the time order operator. In the basis of states one obtains 

<V>ml H (i) I Ipn) = (- 1 )” J„ (L). ( 14 ) 

We will reduce it by the following two steps. Firstly, we can neglect rapidly-oscillating 
terms in e mujt for n > 2. Secondly, we note that | Ji (ic)| | J 0 (:r)| when x is around the 

zeros of |Ji (ic)|, |Ji (x 0 )| = 0. Then when the parameter L ~ x 0 , with x 0 = 2.41, 5.53, 
..., one can drop the term of (-0 m | Hi ( t ) Therefore the approximation and the 

feature of the Bessel function lead to the following approximate expression 


(V’ml Hi ( t ) \4> n ) ^ f (-1)'" Jo ( L ) 6. 


(15) 


i.e., state \if n ) is the simultaneous eigenstate of time-dependent Hamiltonian Hi(t). 
Then in the basis of states {| if n )}, the propagator in Eq. (TT3]) is reduced to 


U Ln {t, 0 ) « Smn exp 


—i 


(i> m \ Hi (t') \f; n ) dt’ 


The corresponding propagator U (t, 0) in Schrodinger picture is in the form 
U mn (■ t , 0 ) ss S mn exp {-i [mix + (-l) m 7 ] t} , 

where 

7=§Jo (L). 


(16) 

( 17 ) 

(18) 


3.2. Dynamics 

To clarify the feature of the dynamics, we consider the time evolution of an arbitrary 
state. We note that any state can be decomposed into even and odd-parity portions, 
which are spanned by the state l^n) with even or odd ?r, respectively. It shows that the 
parity is conserved during the time evolution under an arbitrary (either odd or even) 
H within the approximate framework. We see that the even and odd-parity parts of 
quantum state evolve independently. Furthermore, the appearance of the term H only 
contributes an overall phase on each parts. It is presumable that the profile of the 
evolved state is determined by the factor 7 in a simply manner. It plays a central role 
in the emerging interference phenomena of two subwaves with different parity. 

Actually, for a given initial state 

(19) 

n 

we have 

IV’M) = U (t,0) IV’(O)) 

= e-^ £ e-“/ 2 n \^ 2 n) + £ e-^/^i {^n+i) • ( 20 ) 

n n 

In order to demonstrate the physical picture of state | fjft)), we consider the time 
evolution of the same initial state under the free Hamiltonian Hq, which can be written 
as 

I* (*)> = e- ia “ \i: (0)) = Y. l« ■ 


n 


( 21 ) 
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It is turned out that |^ 0 (£)) exhibits standard Bloch oscillation with period T B = 2tt/oj 
for a local initial state. Furthermore, a straightforward derivation shows that 

Vl> (t)) = cos ( 7 1 ) \ip 0 (t)) - i sin ( 7 1 ) ^0 (t + ^ . ( 22 ) 

It shows that state \if (t)) can be regarded as the superposition of two evolved states 
under the Hamiltonian Hq. In order to get a physical picture of the phenomenon, 
we introduce an operator B w which relates two states \ipo (t)) and \ f>o (t + -rp-)) in the 
following way 

|^o (*)> = '0o (t + ^ . (23) 

Obviously, B n can be expressed in term of H 0 as 

B n = e - iH ° T B/2' (24) 

The physical meaning of operator B n becomes clear if we apply it on the state | k). A 
direct derivation shows that 

B w \k) = exp [iL sin k\ \k + n ). (25) 

It indicates that operator B n is nothing but the 7 r-boost operator that shifts all 
momentum states by n with an extra phase factor. Then if we consider the initial 
state as a Gaussian wave packet (no, ko)) with momentum fc 0 and the center position 
no, approximately, it evolves as 

I- 0 ( 0 ) = cos (yf) e~ lHot \4> (n 0 , ko)) 

— isin (yf) e~ lHot exp (iLsin/c 0 ) 10(n o + n T , k 0 + 7 r)) (26) 

where n n = — Lcosko. Here the Gaussian wave packet has the form 

1 00 2 

\<Hn„,k 0 )) = ^=We-T<”-”.)V‘»”|n). (27) 

v A n=0 

where R is the normalization factor and a determines the half width of the wave packet, 

A = 2\/2 In 2/a in the case of a € 1. We can see that the evolved state represents 
the superposition of Bloch oscillations of two wave packets with a n^-shifted in position 
and a 7 r-shifted in momentum. The amplitude of each oscillation is modulated with 
sinusoidal time dependence. 

Particularly, we see that in the case of the parameters satisfying the equation 

in'l B - (/) e N) (28) 

the wave packet \(f>(n 0 ,k 0 )) is separated as two equal-probability ones, \<f)(n 0 ,k 0 )) and 
its counterpart |0(no + k 0 + n)), at instant t = On the other hand, we see that 
| if (t)) and |ifo (t)) coalesce with each other at instants t — t n , where 
nn , 

tn = —, (n E N) . 

7 


(29) 
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Figure 2. (Color online) Probability transition between two BOs obtained by 
approximate method presented in this paper and numerical simulation. The initial 
wavepacket \i/j (0)) = |0(no,O)) with no = 1.01 x 10 4 and a = 0.1 is involved in the 
Hamiltonian H e g with il = uj and L = 28.89. The black solid line and red dashed line 
are the shapes of Eq. (BUI) . while blue dashed lines represent the Eq. (BT1) . The solid 
blue lines indicate the numerical results for P a ( [t ) and P b (f), which is in accord with 
our prediction that Q can lead to the probability transition in the special intervals. 


To exemplify these features, we give the center positions n a ( t ) , n b ( t ) and 
corresponding probabilities P a (t), P b (t ) of two wave packets as function of time for 
the initial state (0)) = |0(n o ,O)) as the form 

n a {t) — n b + L cos (ut) = n 0 — L sin 2 {pjt/2) (30) 

and 

P a (t) = 1 - P b (t) = cos 2 (yt) (31) 

which is evolved in the Hamiltonian H e g in Eq. (J4j) with constant fh 

However, it is worthy to point out that the above conclusions are obtained within 
the approximate framework, which we will see from the following analysis. One can 
consider the mechanism of the modulated BOs in an alternative way. For a wave 
packet \<fr(no,ko)) at a certain location, the action of Hq is to drive it to evolve as 
a Bloch oscillation, while the action of Hi is to make a transition from \<f> (no, kf)) to its 
counterpart, i.e., 

H 1 \f(n 0 ,k 0 )) ^ \(f)(n 0 ,k 0 + 7r)). (32) 

Meanwhile, the energy difference between two wave packets 
AE = ^ (0 ( n o, k 0 )\H 0 |0 (no, k 0 )) 

- ^ (0 (n 0 , k 0 + tt)\H 0 |0(n o , k 0 + tt)) 

= cos kg (33) 

suppresses this process. It indicates that only under the condition k 0 ~ ±7t/2, such 
a transition has a higher probability to occur due to the fact A E ss 0. We note that 
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the condition k 0 ~ ±7r/2 is satisfied when e~ lHot \(p (n 0 , k 0 )) and e~ lHot \(ft(n 0 , k 0 + tt)) 
overlap with each other in real space. The underlying mechanism of this phenomenon 
is that H i is a local interaction operator. To demonstrate this point, numerical 
simulation is performed by exact diagonalization for the truncated matrix of the 
Hamiltonian H eq rather than The probabilities P a ( t ) and Pb ( t ) are computed by 
\(ijj (0) e lH ote-iHeqt |,0 (o)) |" anc [ |(pf (0) e lH oh+T B /2) e -iH e qt |^ (q ))| 2 respectively. In Fig. 

[ 2 ] we plot the center positions and the probability distributions of two evolved wave 
packets e~ lHot \if (0)) and e~ lH ° ( f + T s/ 2 ) |-0 (o)), obtained by analytical expression in Eqs. 

(1301) and (l3TTh and the numerical simulation. 

It shows that the probability distributions are step-like as functions of time and the 
probability exchange happens when two wave packets meet together. When they meet 
each other for the first time, the splitting of a Gaussian wave packet can be observed. 

It also indicates that our analytical results, Eqs. (T30l) and ()3Tlh is the time-averaged 
approximation for the probability transition. 

A Gaussian wave packet with the center position ft and the half width A = 

2\/2 In 2/a = 23.5, considering the evolution which has been discussed above, can be 
regarded as being local within a spatial extent, D = L + 6 /a, where L = —AgVh/co. 

When ft ~ 10 4 , L ~ 28, D < 100, now we let D = 100; in this case, if we want to satisfy 
| (y/i ~ Vn) j\fn | < 5% for i E D , n ~ 10 3 is mandatory at least. 

In the next section, the validity of the prediction will be investigated in the 
Rabi system by numerical simulation. We will quantitatively evaluate the extent of 
approximation of the above analysis. 


4. Controllable dynamics of cavity-atom system 


Now we apply the obtained results to a concrete case and then demonstrate the dynamic 
property of the system. We investigate the time evolution of the wave packet in the 
single-mode atom-cavity model. As mentioned above, the dynamics of local state in 
large excitation number region is equivalent to that of the Hamiltonian in Eq. (J4j). We 
will show that the modulation of the Bloch oscillations can be controlled by the coupling 
constant and the energy level of the atom, which can be adjusted via external field. 

As an application of the obtained result, let us take a simple case as an example. 
Without loss of generality, we assume that the initial state is in the form 

IV- (0 )) = A±k>^eA(-«q„). (34) 

^ n 

Transforming the basis {|g, n ), |e, ?z)} to {|i) e } and {|i) 0 }, we have 


y(o)> 




/• -\2 
, -5-(*-ra) 


V2R " 

A=e,o 2=0 


\i) 


A > 


(35) 


which corresponds to the superposition of two independent stationary wave packets in 
two chains H^ s with A = e, o, respectively. According to the above analysis, this state 
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evolves to 


IV 1 (*)> = - 7 = e mUlt i & cos ( 7 *) + i^+Tb /2 sin (yt)] \g, n) 

* n 

+ [$ t cos ( 7 1) - id t+ T B /2 sin ( 7 *)] \e, n )} , 


where 


fit = exp 


-ih - y (n - n t y 


h = - y sin ( ut ) > 

- t ■ * 1 
7i t = n — L sin ( — 


(36) 

(37) 

(38) 

(39) 


We focus on the evaluation of the photon number distribution P (n, t) , characterizing 
the dynamics of the state. A straightforward derivation shows that 


P(n,t) = \{g, n \if (t ))| 2 + \{e,n (t ))| 2 


1 r 
R 


l^tl 2 cos 2 ( 7 t) + \d t+ T B / 2 \ 2 sin 2 ( 7 t) 


where \d t \^ has the explicit form 

l^tl 2 = exp [-a 2 (n - n t ) 2 ] , 


(40) 

(41) 


which represents a Gaussian distribution with time-dependent center nt. It is shown 
that the photon number distribution P (n, t) exhibits amplitude modulated BOs. 

We note that P (n,t ) is directly determined by the factor 7 in Eq. (TT 8 lh which 
depends on the ratio gy/K/oj. The numerical method is employed to simulate the 
time evolution process. It is performed by exact diagonalization of the Hamiltonian 
with truncated approximation. We are interested in the suppression and process of the 
probability transition. In the following, we will investigate the cases with two types of 
parameters: i) adjusting the ratio gy/fi/oj for the zero-transition point, and ii) in-phase 
resonant oscillating G (t). 


4-1. Constant G 

According to the approximate analysis, when we take gy/fi/oj satisfying J 0 ( L ) = 0 , the 
transition between two BOs can be frozen. However, in practice a deviation may occur 
with respect to the approximate solution. We compute the time evolution and search 
the zero-transition point for the initial state in the form of Eq. (13411 . We consider four 
typical cases with L = 24.31, 25.73, 27.50, and 28.89. From the plot of the function 
Jo (x) in Fig. (|3]), we can see that two of the four points locate at the vicinity of the 
zeros of the Bessel function, while the others locate at the stationary points. In Fig. 
(141) , the photon number distributions P (n,t) for these parameters are plotted. The 
corresponding analytical results n t (and n t+ T B / 2 ) in Eq. (1391) are also added on the 
plots of P (■ n , t) as a comparison and guide. 
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Figure 3. (Color online) Bessel function Jo (x). Two points a and c with x = 24.31 
and 27.50, are in the vicinity of the zeros. Points b and d with x = 25.73 and 28.89, 
are the stationary points. 




Figure 4. (Color online) Time evolution of the initial state in Eq. (l34l) with the same 
n and a as that in Fig. It is driven by the original Rabi Hamiltonians with f i = uj 
and parameters (see Fig. @) (a) L = 24.31, (b) 25.73, (c) 27.50, and (d) 28.89. The 
photon number distribution P (n, t) is obtained by exact diagonalization. The solid 
and dashed lines are the center positions of two wavepackets from Eq. (15iH) . It shows 
that the probability transition is frozen in the cases of (a) and (c). 
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f.2. In-phase resonant oscillating (t) 


We consider the resonant 11 (i), which have the half period as the Bloch oscillations. To 
clarify the action of If (f), we take a rectangular wave as 

[ ±{nT-T 1 )<\t-<p 0 \<i{nT + T 1 ), 

ft (*) = <"’ (n = 1,2,3...,) (42) 

I 0 otherwise 

with T = Tb, T\ = 0.1 T. Here (po control the phase between 11 (t) and the Bloch 
oscillations. We consider two typical cases with ip 0 = 0 and T/4, which are in- phase 
and out-of-phase with the Bloch oscillations respectively. 

Similarly, we can also take a sinusoidal wave as 


If (t) = 0.5 {1 + cos [2cn (t + ^ 0 )]}a;, 


(43) 


with <po = 0 and T/4, respectively. Here we use a uniform mesh in the time discretization 
for a time-dependent Hamiltonian, i.e., 

r ft 1 


H’ (*)) = exp 


—i 


H (t') dt' 


|^( 0 )). 


(44) 


L Jo J 

The profiles of the evolutions with different types of If (t) are plotted in Fig. [5] The 
numerical result accords with our prediction: in the case of in-phase 11 (t), the wave 
dynamics is similar to the case with constant If, while in the case of out-of-phase If (t), 
is similar to the case with zero If. 


5. Summary 

In summary, we have studied the dynamics of the Rabi Hamiltonian, identifying a 
parameter regime corresponding to amplitude modulated BOs, another type of Bloch- 
Zener oscillations. It is the first time to present a periodic phenomenon for this well- 
known model when the coupling constant goes beyond the RWA. It also reveals a fact 
that a significant effect on the held, Bloch-Zener oscillation of the photon probability 
distribution with a distinct amplitude. It is remarkable that such a macroscopic 
phenomenon is induced by a single atom. Moreover, the probability transition between 
the two BOs can be controlled and suppressed by the ratio g^/ft/to, as well as the in-phase 
resonant oscillating atomic frequency 0 ( t ), leading to multiple zero-transition points. 
The numerical simulation of dynamics in the Rabi model confirms this prediction. 
However, the experimental observation of this prediction requires a coherent time scale 
of g\[fi, which is still a challenge so far. At the same time, with such a high number of 
photons in a quantum resonator, dissipation may have a role, too. 

Acknowledgments 

We acknowledge the support of the National Basic Research Program (973 Program) of 
China under Grant No. 2012CB921900 and CNSF (Grant No. 11374163). 




Amplitude modulated Bloch oscillations of photon probability distribution in cavity-atom systemic 



Figure 5. (Color online) Time evolution of the initial state in Eq. (1M|) with the same 
h and a as that in Fig. @ 4 gy/n/co = 28.89. The fl (t) is taken in the form of Eq. 

with (a) <po = 0 and (b) T/4, and in the form of Eq. (T451) with (c) ipo = 0 and (d) 
T /4. We can see that for both square and sinusoidal waves of S2 ( t ), the profile of the 
evolution in the case of in-phase (out-of-phase) is similar to the case with constant fl 
(zero fl). This accords with our prediction. 
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